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ABSTRACT 


The mean electron density of the solar corona was determined 
hy measuring the dispersion of radiofrequency pulses from pulsar 
UP 0532 during the June 1973 solar occultation. Observations were 
carried out at Arecibo Observatory. Trends continued which were 
noticed in 1971 as solar activity declined from its 1969 and 1970 
levels: 17 r 0 east and west of the sun, 1971 path-integrated density 

was 509s> greater; while 1973 density was approximately 80 % greater 
than the corresponding 1969 and 1970 values; and yet 5 south of 
the sun, 1973 and 1971 integrated densities were only approximately 
one -half their 1969-70 values. A unique solution cannot be found 
by observing the occultation of a single object, but model fitting 
results suggest that the corona continued to become even more con- 
centrated toward the equator in 1973 than in 1971* The equatorial 
number density of electrons at 10 r 0 was inferred to be approximately 
the same as in previous years (7500 cnT-'j, while the number density 
of electrons in all other regions decreased. The decrease of density 
with radius also appeared still less steep in 1973 than in 1971, and 
markedly less than in 1969-70. The coronal multipath scattering 
delay, an index of coronal inhomogeneity, remained near its 1971 
levels, less than one-tenth of its value in 1969-70* 

The best model of the distribution of corona electrons is 
suggested to be one with zero density at the poles, which appears 
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consistent with the results of optical investigations during eclipses 
at a similar phase of the solar cycle. The best model is 
p = 75°° (r/lO r Q ) cos 2 ^ cm ^ , where \ is solar latitude. K- 
corona isophotes and contours of equal path- integrated density are 
presented for several models in «xder that other investigators may 
use their own criteria to Judge which model is the most realistic. 
Ehctensive tables of electron density versus date and position in the 
corona are included for a variety of models. 

A local enhancement of density and a local hole in the corona 
irere observed in 1973. From an analysis of four occultations of 
HP 0532 it is seen that there is no simple relationship between the 
onset of major solar activity and density or scattering enhancements. 
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I. INTRODUCTION 


In recent years the solar corona has "been studied extensively 
by a variety of techniques. Each method has certain advantages and 
limitations. The purposes of this study are to refine the technique 
of determining the coronal electron density by observation of solar 
occultations of pulsars and to provide a means of calibration for 
other techniques. 

A. The Techniques of Measuring Coronal 
Electron Density 

Eclipse Photometry : During a solar eclipse the corona can 

be photographed out to a radial distance of at least 13 solar radii 
[Cox and Liebenberg, 1973 3 • The coronal light consists of two 
components : the K-corona which results from Thomson scattering of 

the sun’s light by free coronal electrons; and the P-corona which 
is caused by the scattering of sunlight by dust in the ecliptic 
plane. In prac:ice, the two components are separated by assuming 
that the emission from the K-corona Is polarized; while the E-corona’s 
light is assumed to be completely unpolarized. 

However, recent calculations show that the F-corona may be 
slightly polarized [Newkirk, 1967]. Beyond 10 solar radii the radiance 
of the K-corona is so small that the inferred electron density would be 
altered drastically if even a minute portion of the P-corna’s light 
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is indeed polarized- Thus the application of this technique is 
presently limited to radial distances less than 10 solar radii. 

Also, even in that region, it is difficult to distinguish coronal 
light from the sky background. 

Radio Source Oc^iltations : The scintillations of continuous 

radio sources as they are occulted by the corona can be used to 
measure the density fluctuations of the corona out to 80 solar radii. 
However, the interpretation of these fluctuations in terms of actual 
coronal density depends upon several uncertain assumptions regarding 
small scale scattering properties of the corona [Newkirk, 1967].. 
such as alignment, scale size, and space-filling factor of the 
density fluctuations. 

Direct Measurements : Spacecraft in heliocentric orbits have 

also carried instruments to directly probe the interplanetary medium 
[Hundhausen, 1972a,]. Presently no in situ measurements have been 
made closer to the sun than Mariner 10 ’s 0.4 a.u. , although the 
Helios spacecraft will also directly measure the medium to within 
0c ‘) a.u. 

Radio Pulse Dispersion : An electromagnetic wave traveling 

through an ionized medium propagates with a group velocity that is a 
function of frequency and of the number density of free electrons. 

The difference in transit times at two frequencies is directly pro- 
portional to the electron number density integrated along the path. 
This technique has been used by spacecraft in heliocentric orbit 
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that receive pulses from the earth at two frequencies [Staff, Stanford 
Center for Eadar Astronomy, I 966 ; Koehler, 1968 ]. Simultaneous radar 
observations of Venus at two frequencies have also been undertaken in 
order to measure the density of the intervening medium [Campbell and 
Muhleman, 1969 ]* Neither of these techniques has yielded information 
on the density of the medium closer than about 0.7 a.u. to the sun. 
Eecently Mariner 10 data have presumably extended the range of dis- 
persion measurements to within 0.4 a.u. Also the Helios spacecraft, 
launched in December 1974, will probe the medium as close as 0.3 a.u. 
from the sun with the dispersion technique. 

Counselman and Shapiro [ 1968 ] suggested that th? dispersion 
technique te used in conjunction with pulsars. It is this method 
that is used in the present work. Pulses from a pulsar are monitored 
as it is occulted by the corona. The relative delays of the received 
pulses at various frequencies can be used to determine the path- 
integrated electron density between the earth and the pulsar. The 
observed electron density is the sum of two components, one attribut- 
able to the corona and the other, to the interstellar medium between 
the observer and the pulsar. 1 The pulsar must be monitored at a wide 
range of solar elongations in order to separate the two components. 

With the exception of the in situ measurements by spacecraft , 
each of the techniques yields only a line -of- sight integral of 

1 The contribution of the earth’s ionosphere to the total 

observed electron density is approximately lO 1 * 1 ' e“ cm“ 2 [Garriott 
et al. , 1955l> which is entirely negligible. 
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electron density. Some assumptions about the geometry of the corona 

I 

are necessary in order to determine the number density at a point in 1 ’ 
the corona. 


B. Earlier Results from the Pulsar 
Dispersion Technique 

The .Jolar oceultation of the Crab Nebula pulsar was observed 
from Arecibo in 1939, 1970 , and 1971* The experiments yielded infor- 
mation on coronal densities in bhe region of 5-20 solar radii. The 
lower limit was fixed by the closest approach of NP 0532 to the sun. 
The upper limit was set by uncertainties in the (time-variable) 
interstellar electron content. The path integrated coronal electron 
density was 3*7 X 10 e" cm at closest approach (5 r ) of the 
earth-pulsar line to the sun in 19^9 and 1970 [Counselman and Rankin, 
19723. For those years, when solar activity was at a peak, the 
corona was assumed to be spherically symmetric. Using such a model 
for the corona, the volume electron density was found to be: 


P « 7000 


r 


10 r 

1 © 



( 1 . 1 ) 


In 1971 [Counselman and Rankin, 19733 the path integrated 
coronal electron density at 5 r Q was only half of its 1969 and 1970 
values. However, four days before and after the closest approach 
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of the pulsar-earth line to the sun, when the pulsar was 17 r Q east 
a.nri west of the sun, respectively, the 1971 integrated density wa s 
50 % greater than the corresponding values measured in 196 9 and 1970* 
The 1971 model-fitting results yielded an equatorial density similar 
to 1969 and 1970 but indicated that the corona was highly concen- 
trated toward the equator and that the logarithmic radial density 
gradient was less steep than in previous years. In 1971 the scat- 
tering due to coronal inhomogeneity also decreased markedly. 

C. The 1973 Experiment 

It was deemed useful to repeat the experiment in 1973 for 
several reasons c First, it is important to continue to monitor the 
changing nature of the corona through varying phases of the solar 
cycle. Second, substantial improvements in the experimental tech- 
niques had been made since 1971 « A newly built digital pulsar anti- 
dispersing device [Boriakoff, 1973 ] had significantly increased the 
S/H available on the Crab Pulsar, and circularly polarized feeds had 
been mounted on the Arecibo radiotelescope in order to eliminate 
the possible timing errors due to the reception of a linearly pol- 
arized pulsar pulse by a linearly polarized antenna whose orientation 
changed through the observing period. Also, two pulsars were dis- 
covered which passed even closer to the sun than DIP 0532 and within 
two weeks of the Crab occultation. Had we been able to observe them, 
the exrbra pulsars would have enabled us to refine greatly the shape 
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of the corona. Finally , R. M. MacQueen asked us to do the experiment 
to assist. in interpretation of coronal observations of the Skylab 
mission. The orbiting laboratory carried a coronagraph which was 
designed to photograph the corona out to 6 solar radii. Because of 
the difficulties in interpretation of photographs of the corona 
mentioned above, it was hoped that our experiment could give a direct 
and unambiguous calibration of the line-of-sight coronal electron 
density recorded on their photographs. 
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II. EXPERIMENTAL DESIGN 


A. Propagation of a Radio Wave in a Cold, 

Tenuous , Homogeneous Plasma 

Let us first consider the properties of a pulse that is 
emitted simultaneously at two frequencies and which propagates 
through a cold, tenuous, homogeneous, and completely ionized medium. 
The quantities derived here form the theoretical framework for the 
experiment . 

Ra nk i n [1970] shows that the difference in arrival times at 
the two frequencies is: 


t 


1 



1 \ 



_1_ 

2c 



path 


ds 


(II. l) 


where and utg are the circular frequencies of the corresponding 
propagating waves and tt is tile electron plasma frequency in the med- 
ium, Rankin showed that terms of higher order which would normally 
appear in Eq. (ll.l) may he neglected when observing NF 0532 in the 
frequency range of this experiment. Equation ( II. l) may he rewritten: 


t 


1 




(II. 2) 
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where D is the dispersion constant, and f.. and f n are the respective 
frequencies . 

The dispersion content is directly proportional to the path- 
integrated electron density, which is called the dispersior. measure 

(dm). 
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(II. 3) 


Thus, In the case of a homogeneous medium, the path-integrated 
electron density can be determined by observation of the difference 
in arrival times of corresponding pulses at two different frequencies 
However, small, irregular density fluctuations in the ionized medium 
between the pulsar and the earth complicate the situation. Their 
effects are discussed in the next section. 

B. Refraction of the Wave Due 
to Density Irregularities 

Eadio scintillation measurements reveal the existence of small 
irregular density fluctuations in the corona. The Irregularities 
randomly scatter the radio waves from the pulsar. Erickson [1964] 
observed the apparent angular broadening of Taurus A as it was occult 
ed by the corona and found the rms scattering angle in radians to be 

cp^ = 1300 f” 2 B~ 2 (II. 4) 


I 
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where R (in solar radii) is the solar impact parameter and f is the 
radio frequency in MHz. The effect of this random refraction is to 
further delay the peak of the pulsar- because of the increased path 
length that must be traversed. The expected scattering delay, then, is 

< T S> = ^ = x 10 "*^ ^ rJ+ sec r ©^ (II. 5) 


where a is the earth-sun distance. 

Thus Eq. (II.2) is no longer strictly correct. The delay at 
each frequency is instead represented by a combination of dispersion 
and scattering effects: 

At “ Df -2 + Sf" 1 * . (II. 6) 

Arrival-time measurements are required at a minimum of three 
frequencies in order to distinguish the dispersion and scattering 
delays. An expression of the form 

*„(») = t o + E 1~ S + Sf" 4 (XI. 6a) 

will be used to model the propagation delays. The observed arrival 
time at a given frequency Is t (f). The three unknown parameters 
t Q , D, and S will be determined. The parameter t is the arriv al 
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time of the pulse if the corona were not present. If the intrinsic 
timing of the pulsar is sufficiently stable, the predicted infinite- 
frequency (corona-free) arrival time may be Included as an additional 
measurements . 


C. Geometry of the Occult at ion 
The path of the sun across the sky in the vicinity of the 
ohree pulsars is shown in Figure 1. 

The pulsar dispersion method, like other line-of-sight techni- 
ques, yields the path -integrated electron density. In order to 
determine the electron density at a point, a model of the solar 
corona must be assumed. In 1969 and 1970 [Counselman and Rankin, 
1972] a spherically symmetric model of the corona of the form 



(H-7) 


was used, p and a were allowed to vary; r was fixed at 10 solar 
o o 

radii. This model must be modified at solar minimum when the corona 
is highly concentrated toward the equator [Newkirk, 1967]. The 1971 
corona was modeled [Counselman and Rankin, 1973] with the addition of 
a term involving the solar latitude k and a "flattening parameter", e: 
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This relation will be used for the 1975 corona. p Q , a, and 
e are allowed to vary. However, as the inclination of the solar 
equator to the ecliptic is only 7.°25, it is impossible in practice 
to distinguish a latitude dependence of density (which would be 
reflected by the value of e) from, a radial distanct dependence (which 
determines a) with a single occultation. 2 The disc- very of pulsars 
JP 05^0 and JP 06ll provided the hope of sampling the electron den** 
sity along three distinct paths across the corona. Had such a 
possibility been realized, the true shape of the corona (and thus e) 
could have been determined much less ambiguously. 


3 The passage of the pulsar south of the solar south pole 
occurs within a day of the closest approach of the pulsar to the sun. 
Thus variations in observed density during the occultation can be 
ascribed equally well to daily changes in the solar latitudes in 
which the (asymmetric) corona is sampled or to the different radial 
distances sampled each day. 
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III. EXPERIMENTAL EQUIPMENT AND TECHNIQUE 

The experimental technique and equipment for the 1973 occulta- 
tion were similar to those used in previous years. The important 
improvements are discussed extensively below, but only an outline of 
basic equipment and technique is given. For a detailed discussion of 
basic equipment and technique see Rankin [1970 3 , Counselman and 
Rankin [1972], and Counselman and Rankin [1973 3 • 

A. The Telescope and Radiometers 

The Arecibo radioteleseope is a 1000-foot diameter, fixed 
spherical reflector with traveling feeds. Any one of approximately 
ten different feed systems can be used by the experimenter. In 
order to differentiate scattering and dispersion this experiment 
requires observation of the pulsar simultaneously at three or more 
frequencies. A new circularly polarized feed system was constructed 
and tested for this experiment in the spring of 1975* This feed 
system replaced the linearly polarized array which had been used 
since the beginning of the Crab surveillance experiment in 1969* 

Thus any systematic errors due to reception of linearly polarized 
pulses by a linearly polarized feed whose orientation changes through- 
out the observations were reduced. 

During the 1973 oecultation experiment the surface of the 
Arecibo reflector was being replaced in order to upgrade its 
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capabilities. While the construction caused a significant increase 
in the power in the side lobes of the telescope as well as minor 
pointing errors, the experiment was not significantly degraded. 

A block diagram of a typical radiometer is shown in Figure 2. 

B. Special Requirements for Observing HP 0552 

The sun passes 5 solar radii north of the Crab Nebula pulsar 
on June 14 or 15. Observations of NP 0532 are complicated by several 
factors. The pulsar is imbedded in the Crab Nebula, a noise source 
whose flux is two orders of magntidue greater than the pulsar's. 

Thus long integration times are required to detect NP 0532. In addi- 
tion, NP 0532 has the shortest period (33 ms) of any pulsar. In 
order to achieve the time resolution required for this experiment, 
an extremely high data rate is required (120,000 samples/second). 

The problems were surmounted through the use of a synchronous sampling 
system. A Doppler -shifted ephemeris for the pulsar's period is used 
to determine a correct sampling rate which assures that consecutive 
pulses are added in phase by the on-line computer. Up to 32,000 
consecutive pulses can be accumulated in this ma nn er. 

C, Interstellar Effects 

The pulsar dispersion technique yields the electron density 
integrated along the path from the pulsar to the earth. It is 
necessary, then, to distinguish the interstellar contribution from 
that due to the corona. Unfortunately, the interstellar dispersion 
and scattering along the path to NP 0532 are time -variable 
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[Counselman and Rankin, 1971; Rankin and Counselman, 1975]. Ko other 
pulsar shows time variability of either scattering or dispersion. ) 
Because of the magnitude and unpredictability of interstellar effects, 
no useful coronal information can be gained beyond approximately 20 
solar radii, where the coronal dispersion and scattering are compar- 
able to the uncertainties in the interstellar baseline. 

In addition to delaying the pulse, the random refractions 
associated with interstellar scattering cause the pulse to become 
significantly broader at low frequencies. (See Figure 3.) Rankin 
and Counselman [1973 1 have studied this phenomenon extensively and 
find that the observed pulse broadening appears to occur in two 
different regions® It Is postulated that one "screen" is the nebula 
itself, while the other consists of the interstellar medium between 
the solar system and the nebula. 

D. Instrumental Effects 

Rankin [1970] found that the pulses are also smeared by the 
receiver system itself. The three contributions are (a) the use of a 
finite bandwidth to observe a pulse sweeping in frequency, (b) the 
rise time of the filter used to limit the IF pas stand, and (c) the 
post-detection filter time constant. The effects can be minimized 
in order to achieve optimum time resolution but the resulting experi- 
mental configuration would require Unacceptably long integration 
times for detection. Such considerations therefore have served only 
as a guide for the actual configuration which is listed in Table 1. 
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E. The Digital Antidispersing Device 

In the spring of 1973 a digital antidisperser he came available 
for pulsar observations at Arecibo Observatory [Boriakoff, 1973]* 

This device electronically adds the proper delays to each of 52 con- 
secutive 20 kHz channels so as to remove dispersion across the pass- 
band. The result is a signal of 640 kHz bandwidth whose characteris- 
tic smearing is only that of a dispersed pulse sweeping through a 
20 KHz filter. A significant increase in S/U results from use of a 
bandwidth much wider than would otherwise be possible. 

F. Allocation of Frequencies 

The antidisperser can be used at any frequency. Considerable 
effort was expended in finding the optimum one. 

The computer can sample each of eight data channels every 64 
jisec. There is considerable freedom in experimental design, with four 
observing frequencies* one antidisperser, and eight data channels 
sampled at 64 jj,sec intervals. The characteristics of a multitude 
of experimental configurations were examined. 

There are several important criteria to be optimized. Cer- 
tainly the most accurate value of the dispersion constant D is highly 
desirable. However, it is also important to minimize the correlation 
between the dispersion and scattering parameters D and S, as well as 
between the dispersion and the infinite -frequency baseline (t Q ). 

It is possible to simulate the various experimental configura- 
tions and determine the expected correlations and standard deviations 


y 
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in D, S, and t without the use of actual data if the expected stan- 
dard deviations of the data are known. Starting witn the expected 
standard deviation in arrival times at eech frequency, the normal 
equations can be inverted to yield the expected standard deviations 
in the parameters D, S, and C, as well as the correlations between 
them. 

Typical standard deviations in arrival times for channels 
without the antidisperser were determined by analysis of NP 0532 
arrival times at each frequency on several days in February 1973* 

The arrival times from l8-minute integrations were referred to a 
common epoch. The rms scatter in arrival times about the mean was 
calculated and called the expected standard deviation. 

The expected improvement in the standard deviation in arrival 
times when the antidisperser is used at any frequency was determined 
through use of the Dicke equation (in which S/N is proportional to 
the square root of the bandwidth), and allowing for smearing in the 
individual 20 kHz channels. 

Table 2 gives the expected standard deviation in arrival time 
for an l8~minute integration at each frequency, with and without the 
antidisperser. (The use of the antidisperser was not considered at 
73.8 MHz because the 20 kHz filters in the antidisperser would indi- 
vidually smear the pulse too much. ) 

Kie expected standard deviations in Table 2 assume a 32 fisec 
sampling interval at 4-30 MHz; 64 p,sec at all other frequencies. (A 
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smaller sampling interval is necessary at 430 MHz due to the intrin- 
sically narrower pulse at that frequency* ) Thus each 430 MHz radio- 
meter actually requires two of the eight channels sampled "by the 
computer at 64 jxsec intervals. 

The results of computer simulations "based upon Table 2 indi- 
cate that the antidisperser is best employed at 196.5 .MHz; with two 
regular channels each at 430.0 and 73*8 MHz, and one regular channel 


at 111.5 MHz. 
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IV. THE CORONAL OCCULIATION 
A. The Pulsars 

The two newly discovered pulsars were not detectable when 
they were occulted by the corona. Thus this experiment was limited 
to an occult at ion of the Crab Pulsar alone. 

B. The Bata 

A typical day’s data consists of from nine to twenty 18-minute 
integrations (’'runs”) at each frequency. The observing configuration 
suggested by the computer simulations was employed until June 13. 

At that time preliminary analysis indicated that little useful in- 
formation was being recorded at 73.8 MHz. Therefore the second 73*8 
MHz channel was replaced by a second 111. 5 MHz channel. The typical 
configuration then consisted of two channels each near 430.0 and 
111.5 MHz; one channel at 73.8 MHz; and the antidisperser at 196-5 
MHz, A regular 196-5 MHz channel occasionally replaced the extra 
111.5 MHz channel in order that the antidispersed 196.5 MHz data 
eouM be directly compared with regular 196.5 MHz data. 

Figure 3 shews actual pulse shapes from l8 -minute integrations 
on July 11 when the sun was remote from the pulsar; and June 15, 
near the time of closest approach of the two. The waveforms are 
significantly noisier on June 15- The dramatic improvement of the 
signal on the antidispersed channel (channel 3) over its conventional 
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196.5 MHz counterpart (channel 7) is obvious on both days. It can 
be seen that the antidisperser provides an extremely well-defined 
pulse shape whose phase is easy to measure. 

C. The Preliminary Analysis 

The first step in the analysis was to determine the arrival 
time of the pulsar’s pulse for each 18-minute run. The experimental 
pulse was cross-correlated with a template pulse whose shape reflects 
both interstellar and instrumental smearing. The maximum of the 
cross correlation function defines the pulse phase and hence the 
pulse arrival time. In addition,, the statistical properties of the 
cross correlation function can be used to estimate the probable error 
in arrival time at 430 MHz. (Experience has shown that there is no 
such relationship at the lower frequencies.) 

Next, a best estimate of the arrival time for the whole day 
was calculated at each frequency in the following manner: The pulse 

arrival time for each run in a day at the given frequency was compared 
to the arrival time predicted by an ephemeris based on the pulsar’s 
past performance. The daily mean of the residuals was then calcu- 
lated. (At 430 MHz the mean was weighted according to the quality of 
the original arrival time estimations, while at the lower frequencies 
all runs were weighted equally. ) The difference between the daily 
mean residual and the residual of the observation made closest to 
transit on that day was added to the measured arrival time of the 
transit r-un. This "phase average" of the daily pulse arrival time 


at each frequency was then weighted according to the scatter of 
residuals within that day. 

Individual runs whose arrival times differed by several stan- 
dard deviations from the average were then carefully examined. Many 
could be reclaimed by re-examining their cross correlation with the 
template pulse and identifying and rejecting noise peaks. At this 
stage the phase of the 73*8 MHz arrivals on each day was found to be 
essentially random. No further analysis of 73*8 MHz observations 
was undertaken. 

The daily phase-averaged arrival times at each frequency were 
next used to make an improved pulsar pulse arrival ephemeris. An 
improved set of daily phase-averaged arrival times was then obtained 
by again comparing each run in a day at a given frequency with the 
predictions of the improved ephemeris. This procedure was iterated 
several times until a highly consistent set of arrival times on each 
day was generated. 

Figures 4, 5* and 6 show the residuals to a quartic fit in 
pulsar phase for 430, 19 &- 5 ) and 111*5 MHz, respectively. Figure 7 
displays the three curves, all plotted with the same seaI Q Data 
from June 8 through 21 were not weighted in the fit. The data points 
represent the "best” arrival time estimate; the error bars represent 
the rms scatter of a day's runs about the mean. The baseline of the 
quartic fit can be taken to represent both the smoothed intrinsic 
timing of the pulsar and timing variations that are caused by 
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changes in interstellar dispersion; any significant deviations from 
it are either due to timing irregularities in the pulsar o: to the 
solar corona. (The characteristic time scale of interstellar disper- 
sion changes is sufficiently long to enable a guartic fit to model 
them accurately. ) It is uncertain whether the broad maximum in 
Figure 4 which lasts from May 26 to July 4 is the result of disper- 
sion delay in the outer corona or of pulsar timing irregularities. 
Because the (corona -free) baseline in Figures 4 through 7 could not 
be established reliably, coronal dispersions could not be determined 
from delays relative to the baseline. Instead, timing differences 
between the various frequencies were used to determine the total 
(interstellar and coronal) dispersion; and the interstellar component 
was removed at a later stage of the analysis. (See below. ) 

Loss of the 75*8 MHz data tended to indicate an increase in 
the strength of the second interstellar scattering screen, which had 
long been quiescent [Rankin and Counselman, 19731. The data at all 
frequencies were therefore cross correlated with templates which 
reflected smearing due to two interstellar screens instead of the 
one screen model which was used in the original analysis. However, 
when the arrival times determin 'd from a double screen model of the 
interstellar medium were used to find daily phase -averaged arrival 
times, th5y were found to be less consistent than those based upon a 
single screen model. Thus the lack of precision in the 73*8 MHz 
measurements can not be explained in terms of increased interstellar 
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scattering. Analysis of recent measurements indicates that the flux 
of WP 0552 at low frequencies has decreased significantly within the 
last few years [Rankin, Payne, and Campbell, 1974]. This is a likely 
explanation of our difficulties at 75*8 MHz. 

D. Path-Integrated Electron Density, 
and Scattering Delay 

A linear weighted least squares fitting technique was used to 
fit the daily average arrival times at 430.0, 196 . 5 ? and 111.5 MHz 
to a model of the form of Eq. (ll.6a). Table 3 and Figures 8 and 9 
display the results. The total dispersion parameter, D, is plotted 
against time in Figure 8. (Recall that 1 sec MHz corresponds to 

]c _ n 

0.743 XIO'e cm .) At this stage it is evident that it is diffi- 
cult to separate the interstellar dispersion from the component of 
dispersion that is attributable to the corona. In Figure 9 the 
scattering parameter S is shown. The only significant scattering 
enhancements occur on June 12 and 14. 

E. The Three-Dimensional Electron Density 
Distribution of the Corona 

1. The Coronal Model 

In the previous section the technique for determining the 
total path-integrated electron density (D) ms described. The para- 
meter D requires no assumptions about the distribution of the corona; 
in this sense it may be considered an observable quantity. In order 
to derive any information on the density at a point in the corona, 


26 


however, one requires some coronal model over which to deconvolve the 
path-integrated variables. The model used in this work, Eq. (II.8), 
includes a dependence upon solar latitude \ and flattening parameter 
e. The exact form of the models dependence upon these two para- 
meters is rather arbitrary. The importance of Eq. (II.8) lies in the 
fact that e and \ allow for any asymmetric density distribution as 
is expected near the time of solar cycle minimum . As Counselman and 
Rankin [1973] point out, such a model should not be taken too liter- 
ally since it is the result of a set of obse-nrations based upon only 
a limited number of rays directed through the corona. While our 
model may not match the details of the actual coronal electron den- 
sity distribution, we believe that its general form permits us to 
perform two important tasks: 

(a) We can subtract the interstellar dispersion from our measure- 
mens in a realistic manner [see (2) and (3) below"!; and 

(b) We can infer at least qualitatively the variation of electron 
density with solar latitude. 

It should be noted that the model does not permit a variation of den- 
sity with solar longitude, and thus implicitly represents an average 
over longitude. 


2. The Interstellar Model 

Equation (II.8) alone is not sufficient to model successfully 
the observed variation of path.-integrated electron density (d) with 
time. The interstellar dispersion was also time -variable. Conse- 
quently the interstellar dispersion was assumed to vary as 
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I>I - B 0 + V* - V + I D 2 (t ’ t o'> 2 (W.1) 

(t was June 15). 
o 

5. Results 

Iterative weighted least squares were used to find simultan- 
eously the three background parameters D q , and D^, as well as 
the coronal parameters p Q and a. All data points listed in Table 3 
and shown in Figure 8 were weighted in the fit. Since our measure- 
ments could not separate a and €, a series of solutions was calcu- 
lated in which e was fixed at various values from 0 (representing a 
spherically symmetric corona) to 1.4 (electron density set to zero 
for all solar latitudes greater than 64,°1 north and south) . The 

results of the various fits are shown in Table 4. In all cases the 

2 

value of the parameter ^ is considerably larger than that expected 

for statistically-independent observation errors. (In 1971 the value 
2 

of x reflected only independent errors.) 

For the 1973 models with e ranging between 0 and 1 . 0 , the 
larger value of \ is ascribable to the significantly improved (~ 0.5x) 
standard deviations in D which are the result of the use of the anti- 
disperser. Thus the 1973 residuals reflect an actual phenomenon 
whose variation Is not being modeled. Nevertheless, no further para- 
meters were added to the model. Figures 10 and 11 show the e = 0.6 
and 1.0 models, respectively, as well as the experimental data. With 
the exception of two large transient variations, the experimental data 
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appear to follow the model adequately. It would be unrealistic to 
attempt to model local density variations in the corona on the 
basis of our line-of -sight measurements. Therefore the present form 
of the model is considered sufficient, and the poor fit Is ascribed 

to (real) local density variations. 

2 

The value of x increases further for models with e greater 
than 1.0. These highly flattened models have no electrons at all at 
polar latitudes. The poorer fit is ascribed to the fact that there 
Is no evidence in the data for the dip in path-integrated electron 
density expected in the extremely flattened models near nid-occultation 
as the pulsar passed near the solar south pole. (See Figure 12, 
which shows the e = 1.4 model as well as the experimental data.) 

An independent means of determining the value of e is obviously 
very important. Unfortunately instrumental difficulties have delayed 
the reduction of Skylab eoronagraph observations. C. F. Keller of 
Los Alamos Scientific Laboratory is also in the process of determining 
K-corona isophotes from photographs taken from an UC-135 aircraft 
during the June 30 eclipse. When these data become available it Is 
expected that the value of e can be determined unambiguously. At 
that time the ambiguous results of Table 4 can be replaced by a best 
coronal model. Until then it Is necessary to use less direct argu- 
ments to determine the value of e. 

First, It is known that the corona becomes highly concentrated 
toward the equator during solar cycle minimum. Photographs from the 
June 30 eclipse as well as Sky lab eoronagraph observations indeed 
depict a highly flattened corona. 
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Second* present models of the corona predict that the solar 
wind velocity is increasing in the region from 5-20 r Q [Hundhausen, 
1972b]. Therefore, as a result of the requirement of conservation 
of mass, values of a less than 2.0 are excluded. This stringent 
criterion rules out models with e less than 0.6, even if a one 
standard deviaticn error in the value of a is permitted. 

Third, when extrapolated to 1 a.u. , all hut the more flattened 
models predict an unrealistically high value for the electron density. 
This criterion must he qualified by noting that our measurements 
were carried out in a regime somewhat different than the interplane- 
tary medium in the vicinity of the earth. 

Fourth, as discussed previously, model fitting results suggest 
that it is unlikely that e is greater than 1.0, although a local 
density enhancement on June 14 could have masked the decrease in polar 
path-integrated density which is predicted hy the extremely flattened 
models hut which was not observed. 

Taken together, the evidence indicates that the corona was 
considerably flattened in June 1973 > with e « 1.0. However, more 
flattened models can not be ruled out. 

Table 5 lists the calculated total dispersion and interstellar 
baseline for e = 0.4, 0.6, 0.8, 1.0, 1.2, and 1.4. The values in 
Table 5 are based upon the fits listed in Table 4. The observed 
coronal electron column densities for a given value of e correspond 
to the difference between the interstellar background on any date 
and the observed total dispersion. Likewise, for a given value of e. 
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the difference between the interstellar background on any date and 
the calculated total dispersion yields a "smoothed" value of coronal 
electron column densities. (The variation of the values of the 
interstellar parameters I> o , and with e in Table 4 makes it 
clear that even coronal column densities require a knowledge of e in 
order to subtract the correct interstellar contribution. Figures 
13 -l8, which are discussed extensively in the next paragraphs, are 
included specifically in order that other investigators can choose 
a best value of e based upon their own criteria.) 

Figures 10-18 portray the best fitting coronal models for 
e = 0.6, 1.0, and 1.4, based upon Table 4. Figures 10-12 show the 
measured coronal dispersion and the best fitting models of dispersion 
versus time for the three values of e. In Figures 13-15 contours of 
equal path-integrated coronal electron density are plotted on the 
basis of the best fitting models for the selected values of e. And 
K-coronal isophotes of the best fitting models for the three values 
of e are shown in Figures 16-18. The isophotes were calculated on 
the assumption that the K-corona results from the Thomson scattering 
of sunlight by electrons distributed according to cur coronal models. 
The inverse -square dependence of the flux of sunlight and an aniso- 
tropic Thomson scattering law were included in the calculations,, 

(Van de Hulst [1950 ] showed that solar limb darkening can be ignored 
safely in the region of interest.) 

The position of the pulsar at Arecibo transit on each day in 
mid- June is also shown in Figures 15-18. It is clear that only a 
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small range of solar impact parameters and position angles has been 
sampled in this experiment. Thus our models are no.t definitive or 
accurate representations of the entire corona. The detailed contours 
shewn for the entire corona in Figures 13-18 are included solely in 
order that other workers may use their own criteria to judge which of 
our models is the most realistic. 

4. A Comparison With Earlier Investigations 
Many groups have done optical studies of the corona during 
solar eclipse. Von KLuber [1958]., Ney et al. [l96l], and Gillett 
et al. [19^43 carried out very careful experiments which should be 
compared with our results, despite the fact that their investigations 
of the K-corona did not extend out as far as 5 r Q , the minimum solar 
impact parameter probed in the present work. 

Von KLuber observed the eclipse of 25 February 1952 from 
Khartoum, when the phase of the solar cyclr was similar to its value 
during the present experiment. On the basis of his intensity and 
polarization measurements, and assuming a spherically symmetric 
corona in his calculations (but not in his results), he presented 
contours of equal electron density out to 3*5 r 0 . The contours 
reveal a very low electron density at the poles, corresponding 
approximately to our e = 0.9 model. 

Key et_ al. observed the eclipse of 2 October 1 959 from two 
stations in French West Africa, when solar activity was near maximum. 
They found that the K-corona appeared slightly concentrated toward 
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the equator even at solar maximum. A model of electron density was 
constructed to fit the observed K- corona in which the density was 
equal to its equatorial value at all latitudes up to some cutoff 
latitude, beyond which the electron density was zero. It was found 
that they could then fit their observations to a corona completely 
free of electrons beyond latitudes of 70 By applying their model 
to von KLuber*s 1952 eclipse data, Hey et_ al. found that no electrons 
were necessary within 25 ° of the poles (i.e. , "beyond 65 ° latitude). 

Gillett et al . , upon observing the eclipse of 20 July 1963 
(also at a solar cycle phase similar to that of the present work), 
also found that a model of electron density with a density cutoff 
near the poles could fit their observed K-corona. 

On the basis of the work of Hey et al. and Gillett et. al. it 
is tempting to select our e = 1.4 model (in which electron density 
is zero beyond 64.°1 latitude) as the best coronal model. However 
our models , in which the electron density is a smoothly-varying func- 
tion of solar latitude, are a compromise between von Kltfber ’ s assump- 
tion of spherical symmetry and the model of Hey et al. and Gillett 
et al. in which the electron density is constant up to some lati - ? .ie 
and then drops to zero. Thus our guess of e « 1.0 does not appear 
inconsistent with von Kluber, Hey et al. , and Gillett et al. 
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F. Local Density and Scattering Enhancements 

The large density enhancement on June 15 (Figures 10, 11, and 
12) may he associated with an importance 2B flare which erupted an 
hour before the start of the Areeibo observations [Carrigan, 19731. 
However, the flare occurred on the solar northern hemisphere (at 
N 15 ° , W 32°), while the pulsar was 5.6 r Q southwest of the sun. 
Surprisingly, no scattering enhancement was observed concurrently. 

The local density decrease observed on June 12 may correspond 
to a hole in the corona, a phenomenon which was occasionally 
observed with the Shy lab coronagraph [MacQueen et al, , 197^3 • How- 
ever, data on that day were taken for only the first half of the 
observing period due to a lightning bolt that struck the control 
room. It is possible that the day's observations are systematically 
biased due to the fact that data were taken only before transit. 

This possibility is considered highly unlikely, however, as all 
known sources of systematic error have been removed, and no such 
bias is observed on days with a full set of data. Although the rms 
scatter of the Jur.e 1? observations is no larger than that of June 
11, the June 12 observations nevertheless represent a less reliable 
estimate of the typical properties of the corona because they do 
not reflect the averaging which occurs over a three-hour observing 
period. Hopefully* Sty lab coronagraph observations will help to 
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determine whether the June 12 point represents a real variation in 
the density. 

The only significant scattering enhancements occurred on June 
12, whose data are slightly suspect, and June 14, the day of closest 
approach, when the earth-pulsar line -of -sight passed almost directly 
below the solar south pole. (See Figure $?. ) At this time coronal 
scattering contributed 3*75 (* 1-0) ms of delay (referenced to a 100 
MHz pulse). 
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V. DISCUSSION AND CONCLUSIONS 


Coronal occult at ions of the Crab Nebula pulsar have now been 
observed over a substantial portion of the present solar cycle. It 
is sueful to compare the results of the four years. 

A- Path- Integrated Electron Density 
Figure 19 is a graphic demonstration of the changes in the 
corona that have accompanied the decline of solar activity from 1970 
to 1973. The peak path-integrated density in 1969-70 was approxi- 
mately twice Its 1971 and 1973 values. But, four days before and 
after conjunction, the 1971 results show 50$ more path-integrated 
density, and the 1973 results show 80 more than 1969-70. These re- 
sults are consistent with a corona whose polar densities are decreas- 
ing as it becomes more concentrated toward the equator. It should be 
emphasized that the curves of Figure 19 are model-dependent. For 
example, if the 1971 value of e was less than 1.0 or the 1973 value 
greater than 1.0, the two curves would tend toward each other at 
Intermediate and large elongations. Or, if less weight is given to 
the 1973 June 15 enhancement in the least squares fits for that year-, 
the peak of the 1973 curve declines. These changes do not alter the 
general conclusion that the corona has redistributed itself toward 
the equator since 1969-70. 
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B, Density Distribution 

The parameters which characterize the overall density distri- 
bution of the corona as inferred from the Crab Pulsar occultation 
experiments in 1969? 1970, 1971, and 1973 are listed in Table 6. The 
equatorial density at 10 r^ has remained relatively constant through 
the four years. The density averaged over all solar latitudes, p Q , 
declined significantly in 1971 and slightly more in 1973- The radial 
density gradient became less steep in 1971 and the trend may have 
continued into 1973* The corona had become highly concentrated 
toward the equator by 1971 and perhaps more so in 1973 , while it 
appeared spherically symmetric in 1969 and 1970* When definitive data 
on the shape of the 1973 corona become available, less ambiguous con- 
clusions can be drawn. However, it does seem reasonable to expect the 
trends which became evident in 1971 to continue into 1973. The Zurich 
smoothed sunspot number, a prime index of solar activity, declined from 
66*7 in June 1971 to 38.6 in June 1973 3 [Solar-Geophysical Data, 1974]. 

C. Radio Scattering 

In 1969 and 1970 the scattering delay in the arrival of a 
pulse referenced to 100 MHz exceeded 1 ms on every day for which the 
impact parameter was less than 12 r^, the reached 34 (± 5*3) ms three 
hours after the eruption of a solar flare. 

3 The 1969 and 1970 occultations , whose results were so similar 
that they were combined, occurred when the Zurich smoothed sunspot 
number was 106.1 and 10. 53* respectively, lending credence to its 
validity as a predictor of overall coronal structure. 
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In 1971 there was no measurable scattering delay except on 
one day, during which a large enhancement of electron density was 
also seen. At that time the scattering delay was 1.4 (± 0.5) ms 
at 100 MHz. In 1973 the scattering delay reached 3*75 (± 1.0) ms 
on the day of closest approach (June 14). No density enhancement 
was concurrently observed. On 12 June 1973 an enhancement of 2.0 ms 
was observed concurrently with a large drop in the density. If the 
coronal void on June 12 is real, our observations of increased scat- 
tering indicate that the matter in the void is highly inhomogeneous. 
This result is consistent with the work of Uoei [ 19731 which, on the 
basis of an analysis of observations of coronal holes, suggests that 
they are the source of high velocity streams of solar wind. Further 
corroboration is provided by Pneuman [1973 1 * who suggests that 
coronal holes are regions of open magnetic field lines which are 
likely to produce large density inhomogeneities. 

Ho other significant scattering delays were observed in 1973. 
It is difficult to make any general conclusions regarding the general 
scattering properties of the corona in 1973 because the data which 
tend to show any scattering delay at all are sparse. Certainly the 
inhomogeneity of the corona has decreased markedly from its 1969 
and 1970 levels. Our data indicate greater inhomogeneity than 
observed in 1971; but such conclusions may be based upon purely 
local coronal phenomena rather than any large-scale properties. 
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D. Local Phenomena 

Density and scattering enhancements were observed on several 

days in 1969 , 1970, 1971, and 1973- (See Table 7- ) The most not- 

l6 -2 

able were an integrated-density enhancement of 10 electrons cm" 
along with an increased scattering delay on 15 June 19^9 associated 

n.(S 

with an active region, an integrated-density enhancement of 10 

-2 

electrons cm on 9 and 10 June 1970 which was identified with a 
narrow coronal streamer at 20 r , an excess scattering delay of '•* 
ms (referenced to 100 MHz) which occurred three hours after the 
eruption of an importance 2 solar flare on 15 June 1970, and an 
integrated density enhancement of 10 electrons cm ~ on l6 June 
1971 tentatively identified with an active region on the back side 
of the sun. Our 1973 results show an integrated density enhancement 
of 10^ electrons cm”^ on June 15 , one hour after the eruption of 
an importance 2 solar flare, and a possible coronal hole and scat- 
tering enhancement on June 12. It is clear that there is no simple 
relationship between density and scattering enhancements. There is 
also no simple relationship between the onset of a major solar flare 
and density or scattering enhancements: The 15 June 1970 flare 

apparently produced a scattering enhancement while the 15 June 1975 
flare apparently produced a density enhancement. 
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Table 1 


Radiometer bandwidths and time constants. 


Frequency 

MHz 

Bandwidth 

KHz 

Post Detection 
Time Constant 
ms 

430.0 

50 

0.1 

196.5 

10 

0.5 

111.5 

5 

1.0 

7:5.8 

5 

1.0 


Table 2 

Expected standard deviations 

in arrival times. 


Expected a 

Expected a 

Frequency 

Without 

With 

MHz 

Antidisperser 

Antidisperser 


lisec 

(isec 

450.0 

120 

35 

196.5 

286 

82 

111.5 

552 

100 

73.8 

3520 



Table $ Total observed dispersion. 

Date = date of observation. 

Solar impact parameter - closest approach of earth-pulsar line 
to center of sun, in solar radii. 

Solar latitude of closest approach = the solar latitude of the 
point of closest approach to the sun of the earth- 
pulsar line. 

Solar equatorial coordinates of earth = X, Y, and Z coordi- 
nates of the earth in solar equatorial coordinates 
(measured in solar radii). 

Observed dispersion = the total (coronal + interstellar) dis- 
persion parameter D. This parameter also corresponds 
to the total path integrated free-electron density from 

the pulsar to the earth, (l sec MHz^ = 0.743 X 10"^ 

- —2 \ 

e cm . ; 
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Table 3 


SOLAR SOLAR SOLAR EQUATORIAL OBSERVED 

IMPACT LATITUDE COORDINATES OP EARTH DISPERSION 

DATE PARAMETER OP CLOSEST ■ - " SEC -MHZ ^ 

R q APPROACH X Y Z 


4/22/73 

169.89 

- 8°2 

159.79 

4/29/73 

152.56 

- 8?5 

176.20 

5/ 9/75 

124.21 

- 9?1 

195.28 

3/12/73 

115.07 

- 9°3 

199.90 

6/ 2/73 

44.54 

-13 ?4 

217.61 

6/ 4/73 

37-42 

-i4?6 

217.92 

6/ 6/73 

30.45 

-l6? 4 

217.99 

6/ 8/73 

23.23 

-19"4 

217.81 

6/ 9/75 

19.73 

-21 °6 

217.63 

6/10/73 

16.27 

-24“8 

217.39 

6/11/75 

12.87 

-29:7 

217.10 

6/12/73 

9.62 

-38?o 

216.74 

6/13/75 

6.78 

-53 f 8 

216,32 

6/14/73 

5.04 

-84 “7 

215.84 

6/15/75 

5.59 

-54?8 

215.30 

6/16/75 

7-93 

-31?3 

214.71 

6/17/75 

10.99 

-19°5 

214.05 

6/18/73 

14.35 

-12“ 9 

213.33 


-147.09 

-18.71 

235734.8 ±3.3 

-127.32 

-16.20 

235740.8 ±2.3 

- 96.ll 

-12.23 

235745.0 ±4.5 

- 86.26 

-10.97 

235741.1 ±1.9 

- 12.80 

- 1.63 

235764.9 ±2.2 

- 5.56 

- O.71 

235769.5 ±1-9 

1.52 

0.x9 

235772.5 + 4.3 

8.86 

1.13 

235771.4 ±3.9 

12.44 

1.58 

235777.8 ± 1.7 

16.01 

2.04 

235779.7 ±1.3 

19.60 

2.49 

235785.0 ± 4.7 

23.18 

2.95 

235776.0 ±3.7 

26.74 

3.4o 

235786.8 ± 2.1 

30.29 

3.85 

235793.6 ±3.1 

33.88 

4.31 

235804.4 ±2.0 

37.39 

4.76 

235793.7 ± 1.7 

40.92 

5.21 

235788.7 ±1.6 

44.48 

5.66 

235788.3 ± 1.6 


Table 3 (cont.) 


SOLAR SOLAR SOLAR EQJJARORIAL OESERVED 

IMEACT LATITUDE COORDINATES OF EARTH DISPERSION 

DATE PARAMETER OF CLOSEST SEC -MHZ 2 

R q APPROACH X Y Z 


6/19/73 

17*76 

- 8°9 

212.56 

6/20/73 

21.26 

- 6°2 

211.73 

6/21/73 

24.76 

- 4°2 

210.84 

6/23/73 

31.82 

- l?7 

208.89 

6/25/73 

38.90 

o;o 

206.70 

6/27/73 

45.94 

+ 1°2 

204.29 

6/29/73 

52.96 

+ 2^0 

201.65 

7/ 1/73 

59*92 

+ 2^6 

198.79 

7/ 4/73 

70.23 

+ 3“4 

194.10 

7/ 7/73 

80. 4o 

+ 3”9 

188.92 

7/ 9/73 

87.04 

+ h °2 

185.22 

7/H/73 

93.6o 

+ 4°4 

LSI. 31 

7/13/73 

100.10 

+ 4°7 

177.18 

7/15/73 

106.44 

+ 4°9 

172.89 


47.97 

6.10 

235785.3 ± 1.8 

51.48 

6.55 

235784.4 ± 1.7 

54.94 

6.99 

235787.3 ± 1.7 

61.85 

7.87 

235785. ± 1.6 

68.70 

8.74 

235790.4 ± 1.7 

75.45 

9.60 

235786. ± 1.7 

82.15 

10.45 

235789.7 ± 1.8 

88.74 

11.29 

235788.6 ±1.9 

98.43 

12.52 

235790.9 ± l.l 

107.90 

13.73 

235793.9 ±1.2 

n4.o3 

14.51 

235795.5 ±1.3 

120.06 

15.27 

235800. ±1.3 

125.98 

16.03 

235801.3 ± 1.8 

131.72 

16.76 

235804.7 ±1.5 


Table 4 


Results of fits to a three -dimensional coronal model, with the flattening parameter fixed- 


e 

(fixed) 

Po 

(10 5 cm" 5 ) 

a 

D o 

2 

s MHz 

Dl 

s MHz /day 

D 2 

(s MHz 2 / day 2 ) 

2 

X 

P EQ 

10 5 cm" 5 

p P 0 LAR 

3 

10 cm 

0 

5.3 

1.4 

0 

0.93 

0. 012 

57 

5.3 

5-3 

0.2 

5.2 

1.5 

2*4 

0.92 

0.011 

56 

5.7 

4,2 

0-4 

5.1 

1.6 

4.6 

0.92 

0.010 

56 

6.1 

3.1 

0.6 

5.1 

1.7 

5-7 

0.91 

0.010 

57 

6.6 

2.0 

0.8 

5-0 

1.9 

8.0 

0.90 

0.008 

58 

7.0 

1.0 

1.0 

5.0 

2.1 

9.6 

0.89 

0. 006 

6l 

7.5 

0.0 

%v 

1* 2^ 

5.0 

2.3 

10.5 

0.88 

0.005 

64 

8.0 

0.0 

1.4 

4-9 

2.4 

11.2 

0.88 

o.oo4 

68 

8.4 

0.0 


e 


= flattening. 


The standard deviations are: 


P o = average coronal density at 10 r . 

a = radial density gradient. 


p = 300 . 

a = 0 . 2 . 



2 

X 


= interstellar dispersion and its first two 

time derivatives (an arbitrary constant 

has been subtracted from D }. 

o 

= sum of squares of post fit residuals, 
divided by standard deviation, with 27 
degrees of freedom. 


D =1.6 for high values of e; 
6 for low values. 

D = 0. 03. 
d 2 = o.oo3. 


p m 

PpOLAR 


= solar equatorial coronal density at 10 r^. 

= solar polar coronal density at 10 r . 

= electron density set to zero for all solar latitudes greater than 70?5 North and South. 


** = electron density set to zero for all latitudes greater than 64.°1 North and South. 


vn 
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Table 5 Calculated interstellar baselines and dispersions, based 

p 

upon the models of Table 4, 235755-0 sec MHz have been 

arbitrarily subtracted from all dispersions and baselines. 

Observed dispersion = the total (coronal + interstellar) 
dispersion parameter, D. 

Calculated dispersion - the total dispersion parameter, D, 

as calculated by the best-fitting models for e = 0.4, 
0.6, 0.8, 1-0, 1.2, and 1.4. 

Interstellar baseline = the interstellar component of the 

total dispersion parameter, as calculated by the best- 
fitting models for e = 0.4, 0.6, 0.8, 1.0, 1.2, and 
1.4. 


Table 5 


o o 
**J Sd 



CBS SHVED 

nATE DX3PE3DI0S 
SEC-ICE 2 

FOR € = 0.4 

CALCULATED IHTERSTEIIAR 
DISEERSIffl BASELINE 

EEC-MHZ 2 SEC-MUZ 2 

FOH e = 0.6 

CALCULATED UTTERS T3LLAR 
CIBFERSIOH BAB ELUTE 

SEC-KHZ 2 SEC-KHZ 2 

FOR e = 0.8 

CALCULATED UiTERETELLAR 
DISESRSIOH babelub 

BEC-tCET EEC -WE 2 

FOH e = 1.0 

CALCULATED UTTEHSTELLAH 
DISFEHSIOH BASSLUiE 

SEC- MHZ 2 SEC-KHZ 2 

FOH 6 = 1.2 

CALCULATED IHTERBTELLAR 
DIEPSR3KU BABELIEE 

SEC-KHZ 2 SEC-HHZ 2 

FOH c 0 1.4 

CALCULATED DTE5RSTE1IAH 
DISEEBSIOH BAEELUIE 

SEC-MKZ 2 BEC-MHS 2 

4/22/73 -0.2 ±3.3 

-0.2 

-2*6 

-0.3 

-2.2 

-0.9 

- 1.9 

-1.4 

-2.0 

-1.8 

-2.2 

-2,0 

-2.3 

4/29/73 5.S ± 2.3 

3.0 

0.3 

3.0 

0.9 

2.9 

1.7 

2.6 

2.1 

2.7 

2.2 

2.7 

2.3 

5/ 9/73 10.0 ± 4,5 

8.6 

5.3 

6-7 

6.1 

9-1 

7.5 

9-5 

—5 

9.7 

9.0 

9-9 

9.3 

5/15/73 6.1 * 1.9 

10.6 

6.9 

10.7 

7.8 

11.1 

9.4 

11.6 

10.5 

11.9 

11.1 

12,1 

11.5 

6/ 2/73 29.9 i 2.2 

29.2 

21.2 

£9.2 

22.3 

29-5 

24.5 

£9.8 

26.2 

30.0 

27.1 

50.2 

27.8 

6/ V73 3*4.5 ± 1.9 

51.9 

22.8 

31-9 

25-9 

32.0 

26 .? 

52.2 

2 ?. 9 

32.5 

28.8 

32.5 

29-5 

6/ 6/73 37.5 ± 4.3 

34-9 

24.5 

34.9 

25.5 

34.9 

27.8 

35-0 

29.5 

35.0 

30.4 

35.1 

31.1 

6/ 6/73 56.'. ±3-9 

38.9 

26.2 

38.7 

27.2 

38.6 

29.5 

3B.5 

31.2 

38.5 

52.2 

5S.5 

32.8 

6/ 9/75 4:-. 6 ± 1.7 

41.2 

27.0 

4i.l 

28.1 

4o.9 

30.3 

4o.S 

32.1 

40.7 

53.0 

40.6 

?5.7 

6'1j0/73 44.7 * 1.5 

44.0 

27.9 

43.9 

2 S .9 

43.7 

31.2 

43,5 

32.9 

43.4 

53-8 

43.3 

34.5 

6/il/ri 50. i 4.7 

47-4 

20.7 

47.4 

29.8 

47,2 

32.1 

47.0 

35.8 

47.0 

34.7 

46.9 

35.4 

6/12/73 Un.a * 3.7 

51.8 

29.6 

51.7 

30.7 

51.6 

33.0 

51.6 

34.6 

51-5 

35-6 

51-5 

36.3 

6/15/73 f’.i. ±2.1 

56.9 

30.5 

56 . 4 

31-6 

56.5 

35.8 

56.1 

55.5 

55.7 

36.4 

55.1 

37.1 

6/14/7 1 31 .,' ± 3.1 

61,0 

51.4 

59.7 

32.5 

58.6 

54.7 

57.2 

36.4 

56.0 

57.3 

55.3 

38.0 

6/15/73 63-4 ± 2.0 

62.3 

32.3 

62.1 

53.4 

62.3 

55.6 

62.5 

37.3 

62-3 

38.2 

61.9 

38.9 

6/16/73 r 8.7 * 1.7 

59-1 

33 . a 

59 . c 

34.3 

60.2 

36,5 

60. 8 

38.2 

61.3 

39-1 

61.8 

39.7 

6/17/73 53.7 ± 1.6 

55.6 

34,2 

56.0 

35.2 

56-1 

3V.4 

56.3 

39-0 

56,4 

4o.o 

56 .6 

40.6 

6/16/73 - -.3 ±1-6 

53.2 

35.1 

53.3 

36.2 

53-3 

38.4 

53.1 

40.0 

53.1 

40.9 

55.1 

41-5 

6/19/73 \3 ±1.8 

51.7 

36.1 

51. B 

37.1 

51-6 

39-3 

51.4 

4o-9 

51-2 

41.8 

51.2 

42.4 

6/SO/75 4 S .4 * 1.7 

50.8 

57.0 

50.9 

38.1 

50.6 

40.2 

50.4 

4l.8 

50.3 

42-7 

50.2 

43-3 

6/21/73 52-3 ±1.7 

50.4 

38.0 

50.4 

39.0 

50.2 

41.2 

50.0 

42.7 

49-9 

43.6 

4g.8 

44.2 

6/23/73 50-0 1 1.6 

50.3 

40.0 

50.2 

4l.C 

50.1 

43.1 

50.0 

44.6 

50.0 

45.4 

49.9 

46.0 

6/25/73 55.4 1 1.7 

50.9 

42.0 

50.8 

43.0 

50, B 

45.0 

50.8 

5 

50.8 

47.3 

50.8 

47.8 

6/27/73 51.0 i 1.7 

51.9 

44.0 

51.8 

45.0 

51-9 

47.0 

51.9 

48.4 

52,0 

49.1 

52.0 

49.7 

6/29/73 54.7 * 1.8 

53.1 

46.1 

53-1 

47-1 

53.2 

49 .O 

53.3 

50.5 

53-4 

51.0 

53.5 

51.5 

7/ 1/7? 53.6 ±1-9 

5^.6 

48.5 

54.6 

49.2 

54.7 

51-0 

54.9 

52-3 

54.9 

53.0 

55.0 

53.4 

7/ 4/73 55.9 ± 1.1 

57.1 

51.5 

57.1 

52.5 

57.3 

54.1 

57.4 

55.3 

57.5 

55.9 

57.5 

56.3 

7/ 7/73 50. Q ±1.2 

59.9 

55.0 

59-9 

55-8 

60.0 

57.3 

60.1 

53.3 

60.2 

58.8 

C .,2 

59-2 

7/ 9/73 6 O .5 ± 1.3 

61.9 

57-2 

61.9 

58.0 

61.9 

59.5 

62.0 

60.4 

62.0 

60.8 

62.0 

6l.l 

7/11/73 65.0 ±1.3 

63.9 

59.5 

63.9 

60.4 

65.9 

6l.6 

63.9 

62,4 

63.9 

62.8 

63.9 

63.1 

7/13/73 66.3 * 1.8 

66.1 

61,9 

66.1 

62.7 

66.0 

63-9 

65-9 

64.6 

65 + 6 

64.9 

65 + 8 

65-1 

7/15/73 69." - - .5 

66.2 

64.3 

68.2 

65-1 

68.1 

66.1 

67-9 

66.7 

67.8 

66.9 

67-7 

67.0 


-p- 

— o 
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Table 6 

A summary of the parameters which characterized 
the overall density distribution of the 
corona in 1969 , 1970, 1971, and 1973. 


Year 

P EQ 

(lo3 em“3) 

Po 

(103 cm“3) 

or 

€ 

1969-70 

7.0 

7.0 

2.9 

~ 0 

1971 

6. 2 -, 8 . 4 

5 . 2-5. 6 

£ 2.0 

£ 0.4 

1973 

6. 6 - 7 . 5 

5. 0 - 5 . 3 

£>2.0 

0.6 


Table 7 

A summary of notable local variations in 
integrated density and scattering. 



Year 

Density 

Variation 

16 - -2 
10 e cm 

Scattering 

Enhancement 

Associated 

Phenomenon 

1969 

June 15 

+ 1 

Yes 

Active region 

1970 

June 9 “10 

+ 1 

Wo 

Coronal streamer 
at 20 r Q 

1970 

June 15 

-- 

Yes 

Solar flare 

1971 

June l 6 

+ 1 

Yes 

Active region (?) 

1973 

June 12 

- 1 

Yes 

Coronal void (?) 

1973 

June 15 

+ 1 

Wo 

Solar flare (?) 
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Figure 1 The path of the sun across the sky in the vicinity of the 
three pulsars. 

Figure 2 A block diagram of a typical radiometer [from Rankin, 1970]. 

Figure 5 Observed pulse shapes based on l8-minute integrations. 
Channel 5 is the antidispersed channel. 

Figure 4 Residuals to a quartic fit in pulse phase at 430.0 MHz. 

The error bars represent the ims deviation of the day's runs 
frcm the mean. (The data for 8-21 June were not weighted in 
the fit.) 

Figure 5 Residuals to a quartic fit in pulse phase at 196.5 MHz* 

The error bars represent the rms deviation of the day's runs 
from the mean. (The data for 8-21 June were not weighted in 
the fit.) 

Figure 6 Residuals to a quartic fit in pulse phase at 111. 5 MHz. 

The eror bars represent the rms deviation of the day's runs 
from the mean. (The data for 8-21 June were not weighted in 
the fit.) 

Figure 7 Residuals to a quartic fit in pulse phase at 430.0, 196.5* 
and 111.5 MHz, all plotted with the same scale. The large 
residuals in mid-June are attributed to the solar corona. 

(The data for 8-21 June were not weighted in the fit.) 
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Figure 8 The total (coronal + interstellar) observed dispersion, D. 
(1 sec MHz 2 = 0.743 X 10 15 e" cm" 2 .) 

Figure 9 The scattering parameter, S. 

Figure 10 The measured coronal dispersion and the best model of 
dispersion versus time for e = 0.6. 

Figure 11 The measured coronal dispersion and the best model of 
dispersion versus time for e = 1.0. 

Figure 12 The measured coronal dispersion and the best model of 
dispersion versus time for e = 1.4. 

Figure 13 Contours of equal path- integrated coronal electron den- 
sity, predicted by the best fitting model for c = 0.6, as 
viewed from the earth in mid- June. The position of the pulsar 
behind the corona is shown at transit for each date in mid- 
June. 

Figure l4 Contours of equal path-integrated coronal electron den- 
sity, predicted by the best fitting model for e = 1.0, as 
viewed from the earth in mid- June. The position of the pulsar 
behind the corona is shown at transit for each date in mid- 
June. 

Figure lg Contours of equal path-integrated coronal electron den- 
sity, predicted by the best-fitting model for e = 1.4, as 
viewed from the earth in mid-June. The position of the pulsar 
behind the corona is shown at transit for each date in mid- 
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June 
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Figure 16 K-coronal intensity isophotes, predicted "by the best 

fitting model for e = 0.6, as viewed from the earth in mid- 
June. Successive controus differ in intensity by a factor of 
5. The position of the pulsar behind the corona is shown at 
transit for eacn date in mid- June. 

Figure 17 K-coronal intensity isophotes, predicted by the best 

fitting model for e = 1.0, as viewed from the earth in mid- 
June. Successive contours differ in intensity by a factor of 
5. The position of the pulsar behind the corona is shown at 
transit for each date in mid-June. 

Figure 18 K-coronal intensity isophotes, predicted by the best 

fitting model for e = 1.4, as viewed from the earth in mid- 
June. Successive contours differ in intensity by a factor 
of 5* The position of the pulsar behind the corona is shown 
at transit for each date in mid- June. 

Figure 19 A comparison of certain best -fitting models of dispersion 
versus time in 1969-1970 (spherically symmetric); 1971 (no 
electrons above the poles); and 1975 (also no electrons above 
the poles). The curves for each successive year have been 
skewed by a quarter day in order that the sun-earth-pulsar 
geometry (rather than the date) is aligned in all three. (The 
1969-1970 and 1971 curves are from Counselman and Rankin [1972] 
and Counselman and Rankin [1973], respectively.) 
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